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We analyse the momentum dependence of static non-uniform susceptibilities of layered local- 
moment systems below Curie (Neel) temperature within the 1/S expansion, the renormalization- 
group approach, and first order of 1 /TV expansion. We argue that the previously known results of the 
spin-wave theory and renormalization-group approach for the transverse spin susceptibility acquire 
strong corrections already at sufficiently low temperatures, which appear due to the interaction of the 
incomping magnon having momentum q with magnons with momenta k < q. Such corrections can 
not be treated in the standard renormalization-group approach, but can be described by both, 1/S 
and 1/N expansions. The results of these expansions can be successfully extrapolated to T — Tm, 
yielding the correct weight of static spin fluctuations, determined by the 0(3) symmetry. For the 
longitudinal susceptibility, the summation of leading terms of 1/S expansion within the parquet 
approach allows to fulfill the sumrule for the weights of transverse and longitudinal fluctuations 
in a broad temperature region below Tm outside the critical regime. We also discuss the effect of 
longitudinal spin fluctuations on the (sublattice) magnetization of layered systems. 



Layered local-moment systems (e.g. layered 
perovskitespj R2MO4 (R is some element, M is 
the transition metal) and undoped parent compounds 
for high-Tc superconductors) are distinctly different 
from the cubic magnets because of the reduced value 
of magnetic transition (Curie or Neel) temperature, 
since it is determined in these systems mainly by weak 
interlayer exchange or anisotropy. Relatively small 
values of magnetic transition temperatures allow for 
both theoretical and experimental study of the evolution 
of magnetic properties in layered systems in the whole 
temperature range T < T/v and some temperature region 
above Tjy- 

At low temperatures magnetic excitations in layered 
(as well, as cubic systems) are described by the spin wave 
theory, considering periodic twists of spins with respect 
to the ordered state. This theory however does not de- 
scribe correctly the thermodynamic properties of layered 
systems in a broad temperature range, since it is applica- 
ble only at the temperatures, which are much lower than 
the magnetic transition temperature. In particular, the 
transition temperatures, predicted by this theory appear 
too large and the critical exponents are not described 
correctly (see, e.g. discussion in Ref. 0). 

This situation is reminiscent of weak itinerant 
magnets [3|, where the Stoner mean- field theory does not 
reproduce correctly thermodynamic properties. For itin- 
erant systems, the shortcomings of Stoner theory had led 
to the formulation of spin-fluctuation theory by Murata 
and Doniach[4(, Dzyaloshinskii and Kondratenko[5], and 
Moriya[|j], who considered the effect of collective excita- 
tions (paramagnons). Recently, an improvement of the 
spin-wave treatment of layered Heisenberg magnets by 
considering the effect of spin fluctuations within the RPA 
(ladder)-type analysis was proposed in Ref. Q. This 
analysis points to a similarity between the two classes of 
the systems, since it shows, that (static) longitudinal spin 
fluctuations in layered Heisenberg magnets, described by 



RPA diagrams, play the role, analogous to the param- 
agnons in weak itinerant magnets. The corrections to the 
sublattice magnetization due to critical (non-spin-wave) 
fluctuations appear however also to be important at not 
too low temperatureSj_as can be explicitly shown within 
the 1/N expansion 

These results motivate to study the effect of different 
type of excitations on the momentum g-dependence of 
the static transverse v H (q, 0) and longitudinal x zz (q, 0) 
spin susceptibility in layered systems (second argument 
corresponds to vanishing frequency). In particular, 
while the momentum dependence of x" 1 (q, iui n ) is ex- 
pected to be dominated by one, three, and higher- 
magnon processes, x zz (q, iui n ) is determined by two- 
magnon processes already in the lowest order of 1/S 
expansion. These two susceptibilities are not however 
fully independent, since they must fulfill the sumrule 
E q , m Jx zz ( ( l. Jw n)+X +_ (q,^)] = S(S+1), which con- 
strains relative magnitude of longitudinal and transverse 
spin fluctuations. In other words, the longitudinal spin 
fluctuations should be also important for the transverse 
susceptibility. 

The spin susceptibility of two-dimensional systems 
was analyzed previously within the Schwinger bo- 
son and modified spin- wave a ppr oaches [9Mlll]. the 
renormalization-group approach jl2l . [l3j and first-order 
1/N expansion 14]. As we discuss in the present paper, 
in the symmetric phase the three former approaches over- 
estimate the spectral weight by a factor 3/2 (see also 
Refs. [j| [l(], Gil), which is important for the fulfill- 
ment of the abovediscussed sumrule. The sumrule ap- 
pears to be fulfilled in 1/N expansion, but violated in 
the renormalization-group and spin- wave approaches. 

Extending the results for the non-uniform magnetic 
susceptibilities to the magnetically-ordered phase allows 
to study both, the regime of low temperatures, where 
the spin-wave approach is applicable, and the critical 
regime, where the results of 1/N expansion can be ap- 
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plied. Such an analysis gives also a possibility to investi- 
gate the crossover between the two regimes, which corre- 
sponds physically to changing relevant magnetic fluctua- 
tions from spin- wave to the critical ones. The sumrule for 
spin susceptibilities can serve as a criterion of the validity 
of the results of theoretical approaches. 

In the following we analyze the results of 1/5 expan- 
sion for transverse and longitudinal spin susceptibilities, 
as well 1/N expansion for the transverse spin suscep- 
tibility, which allows us to understand the momentum 
dependence of these quantities in different temperature 
regimes. 



I. THE MODEL AND THE SPIN- WAVE 
APPROACH 

We consider two-dimensional fcrro- and antiferromag- 
nets with the easy-axis anisotropy, described by the 
Hamiltonian 

H = -\ E J v s ^ - ^E 5 ?^ -dJ2(sd 2 , (i) 

ij ij i 

where J iyi+S± = J,J i>i+ s n = aJ/2 , and J i} = 
otherwise, <5_l and <Jii being the vectors, connecting near- 
est neighbor sites in the same plane and in different 
planes; r\ > and D > are the two-site and single- 
site easy-axis anisotropy parameters. 

The 1/5 expansion can be performed using the Dyson- 
Maleev representation 

St = V2Son , S* = S - a\ ai , (2) 
Sr = V2S(a\ - ^a\a\ ai ), 

where a\,ai are the Bose ideal magnon operators. To 
quadratic order we obtain from Eq. (JT|) in the ferromag- 
netic case (J > 0) the Hamiltonian of interacting spin 
waves 

h = h sw + - ¥ J (qi,q2;q3:q4)(^ 1 a^aq 3 aq 4 ) 

qi...q 4 

x <5qi+q 2 ,q 3 +q4i (3) 

where 

H sw = 2J-E q a q a q , (4) 
q 

Eq = S(J Q ~ Jq) + \J\Sf, 

J q = 2J(cosq x + cosq y ) + a J cosq z is the Fourier trans- 
form of the exchange integrals, / = (25 — 1)D/\JS\ + 
(Jo/J)r]S. The interaction in Eq. ([3]) can be treat per- 
turbatively. 

In the case of a two-sublattice antiferromagnet we sep- 
arate the lattice into A and B sublattices. On the sub- 



lattice A we use the representation ([2]), while on the sub- 
lattice B the "conjugate" representation: 

5+ = V2Sb\ , 5? = -5 + b\bi, ieB (5) 

Sr = V25( 6i - ^btbA), 

where b\,bi are the Bose operators. Introducing opera- 
tors Bq 

Oq = (Bq + S q+Q )/2, 

6l q - (B q - B q+Q )/2, (6) 

where Q = 7r, ...) is the wavevector of the antiferro- 
magnetic structure, up to a constant term, we have the 
Hamiltonian of the same form (j4j, but for the opera- 
tors Bq. Note that in this case Eq in Eq. does not 
have a meaning of an excitation spectrum because of non- 
Bose commutation relations [Bq,B^\ = (5 qp + (5 qiP +Q; 
the true excitation spectrum is determined by the diago- 
nalization of the Hamiltonian with respect to the states 
with the momenta q and q + Q; in the present paper, 
however, we consider only static (classical) contributions 
to susceptibilities and (sublattice) magnetization, there- 
fore commutation relations between operators Bq appear 
to be not important. 

At not too low temperatures (T ^> T q , T q = rJ for fer- 
romagnets and T q = r 1 ! 2 .] for antiferromagnets, where 
r = max(/, a/2)) the result for the sublattice magneti- 
zation in the spin-wave theory reads 

*H*Q>=So-^ln^, (7) 

where So is the ground-state sublattice magnetization, 
qo corresponds to the ultraviolet cutoff parameter in mo- 
mentum space for classical spin fluctuations, which is de- 
termined by temperature for quantum magnets (5 ~ 1): 
qo = [T/(J5)] 1 / 2 in the ferromagnetic quantum case, 
go = T/c in the antiferromagnetic quantum case (c = 
V8\J\S is the spin-wave velocity), and qo = \/32 in the 
classical case (5 ^> 1), see Ref. Q. 

The spin-wave interaction can be treat in the low- 
est order of 1/5 perturbation theory the so-called self- 
consistent spin- wave theory by performing decoupling of 
quartic terms in Eq. ([3]), see, e.g. Ref. @. This the- 
ory yields renormalization of the magnon spectrum: for 
the renormalizcd anisotropy and interlayer coupling pa- 
rameters one finds fx = f(S/S) 2 , ar — aS/S; these 
renormalizations are expected to be qualitatively correct 
outside the critical region (see discussion in Ref. [2]), 
in the latter region the three-dimensional critical fluctu- 
ations due to interlayer coupling or domain wall topo- 
logical contributions due to easy-axis anisotropy are ex- 
pected to be important. For antiferromagnets, one has 
to perform additional renormalization J — » Jj/S in the 
first term of Eq. (j4|), where 7 = 5 + 0.079 is the ground 
state (quantum) renormalization of the exchange inter- 
action (its temperature renormalization can be neglected 



3 



in the considering cases, in the following we also assume 
7 = 5 for ferromagnetic case). In case of consider- 
ing the renormalization of the spectrum within the self- 
consistent spin- wave theory, the remaining interaction in 
Eq. Q should be considered as one-particle irreducible, 
to avoid double counting. 



II. TRANSVERSE SPIN SUSCEPTIBILITY 

A. Spin-wave theory 

We consider first the momentum dependence of 
the static transverse spin susceptibility x H (q> 0) = 
J^ T dr(S^ (t)S~ (0)). Using the representations ^ and 
([5]) and decoupling quartic terms within the spin-wave 
approach, we obtain for both, ferro- and antiferromag- 
nets in the small q limit 



Xsw 



(q,o) 



2S 
\J\iq 2 



(8) 



(for antiferromagnetic case the shift q — > q + Q is to be 
performed) . Note that the result (|5J| can be also obtained 
within the modified spin- wave theory [TTj . The result 
([5]) implies that the static transverse spin susceptibility 
in the spin-wave approach would vanish at the magnetic 
phase transition temperature, which contradicts the re- 
quirement x H (<J) 0) — 2x zz (q, 0) at T = Tm and finite- 
ness of the sum of averaged squares of spin components, 
see discussion below, in Sect. IIIB. 

To overcome the latter drawback of the spin-wave ap- 
proach, we consider the lowest-order correction to the 
result ||5J) due to the magnon interaction 



Xc-sw 



(q,o) 

T 2 1 



XsV(q>o) (9) 
v(q,qi;q2,q + qi -02) 



(JlY Q 2 (q 2 + r)[qt + r)((q + qi - q 2 ) 2 + r) ' 

where 'c-sw' stands for the spin-wave result, corrected 
with account of the spin-wave interaction. For the case 
of vanishing interlayer coupling we have 

V 3 (qi) q2; q3> li) = ^q 3 + J<a — ^qi-q 3 — ^qi-q4 

~ -2|J|(qiqa + /). (10) 



The calculation of the integral yields 

X(tsw(q 5 o) = 



s + J-f^h,* 



(ii) 



\J\W 



So \2tt\J\ 7 /1/2 



This result keeps its form also in the presence of the inter- 
layer coupling with the replacement / — > r. We would like 
to stress that the term proportional to ln 2 ^// 1 / 2 ) comes 
from the integration over momenta f 1 / 2 < k < q, and 
therefore inaccessible for the standard renormalization- 
group techniques [l2l [l3| . which consider only scales, 



which are larger than all the infrared cutoffs, i.e. k > 
maxfe/ 1 / 2 ). 

At the same time, as we argue below, the second term 
in Eq. (fTTj) appears to be crucially important for fulfill- 
ment of the sumrule already by the lowest-order pertur- 
bative correction. Indeed, considering the sum of squares 
of the transverse spin components for the result (|lll) we 
obtain 



<(£? ) 2 )stat + ((^) 2 )stat =T £ XctswM) 



r 1 ' 2 <q<q 



2S(S -S) + - W (S ~S) 3 . 
000 



(12) 



The result (TT21 describes correctly not only the low- 
temperature behavior of the weight of transverse fluc- 
tuations, but also the limit T — > Tm = Tc(Tjv). Indeed, 
in this limit we have S = and Eq. HHJ) yields (2/3)5q 
which, as we see in the following, is the weight, required 
by the 0(3) symmetry. 



B. 1/N expansion and comparison of different 
approaches 

To treat the momentum-dependent transverse suscep- 
tibility beyond lowest order of l/S'-perturbation theory, 
we perform 1/N expansion in a way, which is similar to 
that described in Ref. [lH for 2D antiferromagnets (see 
Appendix A). The result of the corresponding expansion 
outside the critical regime is 



>-l/JV,: 



\J\jq 2 



S- 



(271-7! J I ^ r^ 5 ) 



s- 



t m q 

7T7|J| r l/2 



(13) 

It has the same structure, as the result of the l/S'- 
perturbation theory (jlip . except for the denominator in 
the second term, which accounts for the effect of longitu- 
dinal spin fluctuations, and can be understood in terms 
of the vertex corrections, discussed below in Sect. IIIA. 
On the other hand, at T — > Tm we obtain 



^l/N,T=T, 



aT In- 9 



\J\jq 2 V2tt|J|7 rV2 



, q > r 



1/2 



(14) 

with a = 2/3, which is similar to the two-dimensional 
case[Hj]. The results (p2J) at T -> T M and O differ 
only by a coefficient in front of logarithm: Eq. (|13l) cor- 
responds to a — 1/2. This difference appears due to the 
fact, that apart from the spin- wave and longitudinal ex- 
citations, the result of the 1/N expansion in the critical 
regime (fill accounts for non-spin-wave (critical) fluctu- 
ations, which change the coefficient in front of the loga- 
rithm to a = 2/3. 

From the result (|13|) we find the weight of the trans- 
verse spin fluctuations in 1/N expansion outside the crit- 
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FIG. 1: (Color online) The weight of the transverse spin 
fluctuations vs. (sublattice) magnetisation in different ap- 
proaches (from top to bottom): renormalization-group ap- 
proach, Eq. (|17[) . corrected spin- wave theory, Eq. (|11[) . 1/N- 
expansion outside the critical regime, Eq. (|13|) . and standard 
spin- wave theory, Eq. @. The semi-circle at the left axis 
shows the value of the weight a — 2/3, required by 0(3) sym- 
metry. 

ical regime 

<(S?) 2 W + ((^) 2 )stat = r^ x &(q,o) 

q 

= (5o/2)[CT + l-2 < T 2 + cr 2 arctanh(l-cr)], (15) 
where a = S/S , while from Eq. ((HI) we find at T -> T N 

((^) 2 ) s tat + <(S?) 2 W = T^ X +- (q,0) = aS*. (16) 

q 

The weight of the transverse spin fluctuations in the 
discussed approaches is shown on Fig. 1. One can see 
that the weights of the corrected spin-wave theory and 
1 /N expansion outside the critical regime are numerically 
close to each other. Both these theories yield also rea- 
sonable extrapolation of the spectral weight to the limit 
T — > T\i (a — > 0); the momentum dependence of the 
result of 1/N expansion, Eq. ([P31 is expected to be more 
correct in the intermediate temperature region, since it 
qualitatively reproduces also the limiting momentum de- 
pendence, Eq. ([HI) at T -> T M . 

The results (fTTI) and (fT4"f can be compared to the re- 
sult of the renormalization-group analysis (l2[ li~3l. Ii~5| . as- 
suming that the logarithm in the numerator of Eq. (jHJ) , 
contained within the spin- wave theory in S, is in fact cut 
at the largest low-energy scale, which yields 

\J\1T \ 2tt| J| 7 q J 

The result (|17p can be put in the symmetric phase 
(5 = 0) into form (|14p with org = 1. Essentially the 
same result was obtained also in the paramagnetic phase 



within the Schwinger boson spin- wave approach [9|, [lOj . 
(Note that the modified spin- wave theory [ll| also yields 
the result (fPTj) . but for the longitudinal spin susceptibil- 
ity; as discussed above the transverse spin susceptibility 
vanishes in the latter approach in the symmetric phase) . 

One can see, that the low-temperature results of sys- 
tematic first-order 1/5* expansion (fTTj) . as well as 1/iV 
expansion, Eq. (I13|) are distinctly different from the re- 
sults of RG or Schwinger boson approaches, as they con- 
tain second power of logarithm instead of the first. This 
leads to a difference of the weight of the transverse spin 
fluctuations of the abovementioned approaches already 
at sufficiently low temperatures (a close to 1), see Fig. 
1. On the other hand, in comparison with the result of 
1/N expansion at T = T N , Eq. jH]), the result ([17]) 
yields different coefficient org = 1 in Eq. (IT41 . which 
yields incorrect weight ((5f ) 2 W + ((Sf ) 2 ) sta t in the RG 
approach at T ->• T M [13]. 

III. LONGITUDINAL SPIN FLUCTUATIONS 

A. Spin-wave and RG analysis of longitudinal spin 
susceptibility 

Let now analyze the longitudinal spin susceptibility 
within the spin-wave and RG approaches. The spin-wave 
theory yields for the spin susceptibility, longitudinal with 
respect to the order parameter, the result, which is ex- 
pressed as a convolution of the two magnon Green func- 
tions. In particular, in the 2D case with the easy-axis 
anisotropy we obtain 

X0iq) = (J^Y E (P 2 + /T)[(p~q) 2 + /T] 

„ f T/[2^(J 7(7 ) 2 ]ln(g 2 // T ), g 2 » f T 

-\ X = T/[47T(J 7 ) 2 / T ], q 2 « h ■ U8J 

Similar result with the replacement fx — > cut can be 
also obtained for weak interlayer coupling. In compar- 
ison to the transverse susceptibility, Eq. ([5]). equation 
(|18p has an additional log-singuliarity in the small mo- 
mentum limit. As discussed in the previous Section, they 
are not sufficient to describe spin fluctuations at not too 
low temperatures. In particular, the result (|18p does not 
account for the interaction of the spin-waves. 

To treat the longitudinal spin fluctuations beyond spin- 
wave theory, the ladder approximation was suggested @. 
This approximation treats contributions to the vertex 
of magnon interaction in a certain channel, correspond- 
ing to the scattering of two magnons with close mo- 
menta ("longitudinal particle- hole" channel). Below we 
argue that this approximation can be rigorously justi- 
fied by considering logarithmically large contributions in 
the parquet-type diagrams, which appear due to specific 
momentum dependence of the susceptibility (|18l) . 

To perform summation of these contributions, we con- 
sider renormalization of the magnon interaction vertex 
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by magnon scattering in three different channels. In the 
lowest order of the perturbation theory we obtain 



*(k, p - q; k - q, p) = - 2 Jk(p q) 

, T ,2V t k ( s ~ q)Mp- q)l 

S 1 

- at j 2 v (ks)[(s + k - p)(p - q)l 

S 

+ 4TJ 2 



(19) 



k(p-q)[s(s-k-p + q)] 



s-^s — k— p+q 

Using the relation £] s ( s a ~ g a )sb/(E s E s _ q ) = 
5 a b/(AiT)hx(qo/q), we find to second order in 1/5 1 



4>(k,p 
1 + 



-q;k- 
T 

2n\J\ r 



q,p) = -2Jk(p 

-Mil- 

9 |k-] 



(20) 



In — 



In- 



9o 



To go beyond second order of the 1 / S'-perturbation the- 
ory, we apply the one-loop renormalization-group ap- 
proach, which is also equivalent to summation of loga- 
rithmic divergencies of the parquet diagrams. To this 
end we replace infrared cutoffs in the particle-hole and 
particle-particle contributions in Eq. ((20)) by the cutoff 
parameter ^i, differentiate r.h.s. of Eq. (|20|) over and 
replace the bare vertices in the right-hand side of equa- 
tion (l20l) by renormalized vertices [we assume that the 
momentum dependence 

$(k, p - q; k q, p) = -2 J(f>k(p - q) 

does not change, as follows from the lowest-order correc- 
tion ([20])]. In this way we obtain 



^ d/i 



T 



2tt|J| 7 2 



(21) 



The equation (|21[) was earlier obtained within the 
Holstein-Primakoff representation of spin operators in 
Ref. 18]. Using of the Dyson-Maleev representation 
makes however transparent that due to the change of sign 
in the particle-particle channel, which occur due to bilin- 
ear structure of the interaction vertex, the contribution 
of the two of the three channels cancel each other. 

For q ~ |k — p| ~ |k + p — q| we stop the flow at 
fi = max(<7, |k — p|, |k + p — q|) to obtain the result 



1 



1 



2tt| J\t 



ln- 



go 



(22) 



<9,|k-p|,|k+p-q|; 



When the momentum transfer in one (some) channel(s) 
is much smaller than in the others, one has to continue 
scaling with the contribution of the channels with smaller 
momentum transfer retained. For q <§C |k — p|, |k + p — q| 
we obtain again the equation (|22|) with the replacement 
max(g, |k — p|, |k + p — q|) — > q. On the other hand, at 
q |k— p|, |k+p — q| the result (|2"2")l remains valid, since 



the contributions of the two remaining channels cancel 
each other. Therefore, for arbitrary q and |k — p| ~ 
|k + p — q| we obtain with the logarithmic accuracy the 
result (cf. Ref. 0) 



$(k,p-q;k-q, p) 



_ |J|k(q-p) 

S/(2 7 ) + <z 2 xo(<z) 



0(\J\f), 



(23) 

where S is given by Eq. ([7]). Thus, as well as in RPA 
for itinerant magnets j3j, the effective interaction is en- 
hanced by fluctuations. 

For the static (staggered) non-uniform longitudinal 
susceptibility (for antiferromagnetic case the shift q — > 
q + Q has to be performed) we obtain Q 



(S/S )xo(q) 



1 - (r/2 7 r|J| 7 5o)ln[ 9o /max(A 1 /2 ;(? )] 

Sxo(q) 



S +(|J| 7 /2)9 2 Xo(9)' 



(24) 



The vanishing of the weight of longitudinal fluctuations 
at T ~ > Tj\/ is the unphysical result of the parquet ap- 
proach, which do not consider critical fluctuations in the 
vicinity of the magnetic transition temperature. For the 
actual layered systems the range a < 0.5 corresponds 
to the critical regime, which occur in the rather narrow 
temperature range near the transition temperature (the 
magnetization drops rather sharply near Tu, see Sect. 
IIIC). The description of this regime requires consider- 
ing critical fluctuations of either d = 3 0(3) or d = 2 
Zi symmetry, depending whether interlayer coupling or 
anisotropy dominates, see Refs. @, i-Q. At the same 
time, the results ([24]) and (f27|) are expected to be valid 
outside the critical regime, i.e. for a > 0.5. 



B. Sum rule and the weight of longitudinal 
fluctuations 



In the following we analyze the fulfillment of the sum 
rule 



tat + ((^) 2 ) stat + ((S 2 ) 2 ) st at (25) 

= 5> + -(q,0)+ X 2Z (q,0)]=^ 
q 

by the considered approaches; 'stat' stands for contri- 
bution of classical spin fluctuations (with the ultraviolet 
cutoff, given by go). This sumrule follows from the stan- 
dard relation (S 2 ) = S(S+1); the difference S(S+l)-~sl 
in the right-hand side appears due to retaining only clas- 
sical term in equation (|25[) . The possibility of considering 
separately classical and quantum contributions is realized 
for layered systems at temperatures T q <C T \ J\S 2 and 
reflects the 'frozenness' of magnetic moments (i.e. sepa- 
ration of the contributions of quantum and classical spin 
fluctuations) in the state with long-range or strong short- 
range order. The sum rule (I25|) can be rigorously proven 
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within the nonlinear-sigma model treatment of the model 
(fTj) by excluding the dynamic spin fluctuations which im- 
plies quantum renormalization of the model parameters 
similarly to the 4 model analysis see also Refs. 

It can be easily verified that the results of the spin- wave 
theory © and (JTHJ) fulfill the sumrule (US]). As discussed 
in previous sections, this theory is however expected to 
be insufficient at not very low temperatures. Fulfillment 
of the sum rule ([21)1 by the other approaches, consid- 
ered in Sect. II, requires the weight of the longitudinal 
fluctuations to be 



[<(S* 



-s 2 ysl 



(26) 



c-sw 
l/N,non-crit 
l/N,T->T M 
RG 




(l-af(l-l(l-a)) 
i[l — a — a 1 arctanh(l — a)] 
1/3 


2 

The substruction of S in the left hand side of Eq. ([26]) 
corresponds extracting the uniform (Bragg) contribution, 
the result therefore reflects the contribution of the lon- 
gitudinal fluctuations. As discussed above, only the cor- 
rected spin- wave theory and the l/N expansion fulfill 
the 0(3) symmetry, which require that at T — > Tm 
the weight of the transverse components should be twice 
bigge than the longitudinal ones. For the scaling analy- 
sis, the whole spectral weight (excluding the Bragg con- 
tribution) is in fact entirely contained in the transverse 
part (see also Ref. |15j). which explicitly violates 0(3) 
symmentry at T = Tn. 

To verify the conformance of the parquet (RG) ap- 
proach to longitudinal fluctuations to the approxima- 
tions, listed in Eq. (|26[) . we perform summation over 
momenta of Eq. (pM|) to obtain: 

T^ x 2Z (q,0) = 2S[S -S + Sln(S/S )) 
q 

f (So-S) 2 (l-(2/3)(So~S)), 

\ 2S S S -» 

Comparing the first row of the second line of Eq. (|2"T)l 
with Eq. (l26l) we find that at low temperatures it 
agrees with the corrected spin- wave analysis of the trans- 
verse part implying that the considered parquet ap- 
proach yields correct asympthotic behavior of the spec- 
tral weight at low temperatures. In fact, the correspond- 
ing range of sublattice magnetizations a > 0.4, where the 
results for the weight of longitudinal fluctuations in the 
parquet analysis, corrected spin- wave theory and l/N 
expansion outside the critical regime agree with each 
other, appears to be rather broad (see Fig. 2); as dis- 
cussed above, it corresponds to the temperatures below 
the (narrow) critical regime near Tm- As discussed in 
Sect. Ha, the result of the low-temperature l/N expan- 
sion (|13p also allows to fulfill approximately the sumrule 
at a > 0.4; while the result (fT4"[) fulfills the sum rule at 
T = T M . 



FIG. 2: (Color online) The weight of the longitudinal fluc- 
tuations vs. (sublattice) magnetisation. Lines from top to 
bottom correspond to: standard spin- wave theory, Eq. (|18|) . 
applying the sumrule ()25[) to the results of the 1/JV-expansion 
outside the critical regime and corrected spin- wave theory (see 
Eq. ([26])), and the weight obtained in the parquet (RG) ap- 
proach, Eq. ((24)) . 



C. The effect on the (sublattice) magnetization 

Finally, we discuss the effect of longitudinal static spin 
fluctuations on the (sublattice) magnetization. It was 
argued in Ref. 0] that the longitudinal spin fluctuations 
yield the subleading logarithmic correction to the spin- 
wave result (f7]) , which improves comparison of the result 
with the experimental data[2lj]. Such a correction can be 
also naturally obtained in the two-loop renormalization- 
group approach Q, which yields 



1 



T 



(27) 



4irp s 
-2 In 



(N - 2) In 



f 



4 

A(T) 



So 




max(cr 1 //32 5 T/(47r / 9 ;s )) 
-2(1 - a 1/p2 ) + F (T/{A^p s (j 1/h ) 



(28) 



where p s is the ground-state spin stiffness (e.g. p s — 
J 7 So in the self-consistent spin- wave theory), A(T) = 
max(ar/2, /t), and F (x) is some non-singular function, 
which accounts for the contribution of higher-loop terms. 
The "critical exponent" fa, which is the limit of the ex- 
ponent fa+e in d = 2 + e dimensions at e — > 0, is given 
by fa = (N - 1)/(2(JV - 2)). At N = 3 we have fa = 1. 
The leading logarithmic term in Eq. (|28j) corresponds in 
this case to the self-consistent spin- wave theory, while the 
subleading logarithmic term coincides with the result of 
the parquet approximation, discussed in Sect. IIIB (cf. 
Ref. Q), and therefore describes the contribution of the 
static longitudinal fluctuations to (sublattice) magneti- 
zation. The nonsingular term —2(1 — a) docs not follow 
directly from the parquet approach and at low tempera- 
tures compensates the second term in the square brack- 
ets. The function F stands for the contribution of non- 
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parquet diagrams, not accounted by renormalization- 
group analysis. 

The result (|28|) can be compared to that of the first- 



order 1/N expansion of the sublattice magnetization, 
which reads 



expa 
[6, 7] 



T 



(N - 2) In 



2T 2 



-2(1 



4irp s 

2n t , ^^^02 



A(a r , f r ) 



B 2 ln\ (29) 



where a r — «t=o and f r = /t=o are the ground-state 
quantum-renormalized parameters of the anisotropy and 
interlayer coupling, B 2 = 3 + f r / yj + 2a r f r , 



4wp s 



(N-2)T 



(30) 



For the cases, when only anisotropy or interlayer coupling 
is present (a r = or f r — 0), the results (|2"5)| and (|2"9)l 
differ by the replacement ct 1 /' 32 — > a 2 in the right-hand 
side, which is related to the property of the first-order 
1/N expansion, that the N — oo value of the exponent 

= 1/2 is corrected only in the leading (proportional to 
N — 2), but not in the subleading (proportional to N°) 
terms in Eq. (f2"9")l . Although for the abovementioned 
cases the result (|28|) can be therefore considered as a 
self-consistent modification of the first-order 1/N result, 
it appears that non-self-consistent equation (|29j) better 
agree with experimental data (see below), which can be 
related to the importance of non-spin-wave fluctuations, 
not captured by the renormalization-group approach. 

Let us also discuss the relation of the results of 1/N 
expansion for sublattice magnetization in the presence of 
interlayer coupling only (f r = 0) to the approach, pro- 
posed recently in Ref. [l9j |. The result of Ref. [l9[ up to 
the logarithmic accuracy reads 



(7=1- 



T 
2irp s 



hi 



2T 2 



c^aa 



1/2 



(31) 



This result can be compared to the results of first- 
order 1/N expansion ([29|l . One can see, that although 
Eq. (|3"TT) contains some of the first-order 1/N corrections, 
it borrows the exponent of magnetization 1/2 from the 
result of the zeroth-order 1 /N expansion in the interme- 
diate temperature regime and neglects subleading terms 
in Eq. (|29p , which account for the longitudinal spin fluc- 
tuations. 

Comparison of different approaches is presented in Fig. 
3. We choose the parameters, which correspond the com- 
pound La 2 Cu0 4 : p s =290K, c =2618K, and a = KT 3 . 
One can see consequent reduction of the Neel tempera- 
ture from spin-wave and self-econsistent spin-wave anal- 
ysis, and further to the renormalization-group approach. 
The results (|29| and (|3T1) appear to be numerically close 
to each other, and therefore Eq. (|3"Tj) provides success- 
ful extrapolation of the sublattice magnetization from its 
low-temperature limit. 




100 200 



300 400 
T, K 



500 600 



FIG. 3: Temperature dependence of sublattice magnetization 
of La2Cu04 in different approaches (from top to bottom): 
spin- wave theory, Eq. (0, self-consistent spin- wave theory, 
Eq. ( [7J with the replacement a —¥ <xt, renormalization- 
group approach, Eq. (|28[) , 1/N expansion, Eq. ((29}, and the 
approach of Ref. [1S|], Eq. (|31[) of the present paper 



IV. CONCLUSION 

In the present paper we have considered momentum 
dependence of static transverse- and longitudinal spin 
susceptibilities, analyzed the fulfillment of the sumrule 
(f2"5f and discussed the effect of longitudinal spin fluctu- 
ations on the (sublattice) magnetization. 

For the transverse susceptibility, we have shown that 
the result of the corrected spin-wave theory, account- 
ing for the spin-wave interaction in the lowest order in 
1/5, yields the correct weight of transverse spin fluctu- 
ations in both, the low-temperature limit ant T — > Tm- 
The correction to the pure spin-wave result is propor- 
tional to ln 2 ((j , /r 1 / 2 ) and comes from the integration over 
the range r 1 / 2 < k < q of the momenta k of virtual 
magnons, interacting with the physical state at momen- 
tum q. This momentum range is inaccessible for the stan- 
dard renormalization-group approach, which treats only 
contribution of momenta k > q. Therefore, standard 
renormalization-group approach by Nelson and Pelkovitz 
12] and Chakraverty, Halperin, and Nelson ap- 
plied to the non-uniform spin susceptibility, has to be 
supplemented by considering the RG flow at the scales 
k < q. This can be done in particular within the func- 
tional renormalization-group analysis [22j , which is leaved 
as a subject for future studies. 

The first-order 1/N expansion allows to determine the 
static transverse spin susceptibility outside the critical 
regime and in the limit of Curie (Neel) temperature. 
The corresponding correction to the spin-wave result 
for the transverse nonuniform spin susceptibility changes 
from ln 2 (g/r 1 / 2 ) at low temperatures to ln^/r 1 / 2 ) near 
T = T M - The result of 1/N expansion at T = T M dif- 
fers by a factor 2/3 from the results of Schwinger-boson 
and the renormalization-group analysis, which happens 
because the latter approaches assume spin-wave picture 
of the excitation spectrum instead of a critical one. Fur- 
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ther analysis of the crossover between the spin- wave and 
critical regime seems to be of the certain interest. 

The applicability of the ladder (parquet) approxima- 
tion for the longitudinal spin susceptibility was analyzed. 
It was shown that the results of these approximations al- 
low to fulfill the sumrule at the temperatures outside the 
critical regime, where the relative (sublattice) magneti- 
zation a > 0.4. Extending the results of 1/N analysis 
to the order 1/N 2 , which would make possible to de- 
scribe both, longitudinal and transverse spin susceptibil- 
ities within this approach with higher accuracy, are of the 
certain interest. The extension of the presented analysis 
to calculate the dynamic magnetic susceptibilities on the 
basis of the discussed approaches would allow to study 
the difference of the spin- wave and critical dynamics. 
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of longitudinal spin fluctuations to the transverse spin 
susceptibility, as n _1 (<7) has a very similar structure to 
the vertex (|23|). renormalized by the longitudinal fluctu- 
ations. 

Evaluation of the integral in (|34|) yields 



k 2 



Hi _L 

27T /V2 



Hk/f 



l/2\ 



Nln(k/f 1 / 2 ) + 2TT<T 2 /{Tg) 

k 2 1 

~ JV Hn(fc//V2) + 2Tra 2 /(Tg) 



)• (35) 



Using the result of the zeroth-order 1 /N expansion 

( j 2 = l-(gT/(27r))ln(q /f 1 / 2 ) (36) 

(which is consistent with the calculating the first-order 
1/JV correction to the self-energy), we represent Eq. (|35|) 
in the form 



Appendix A. 1/N expansion for the transverse spin 
susceptibility 

Following to jf| 0, [2(|, to construct 1/N expansion 
we pass to the continuum classical 0(JV)-model with the 
partition function 



l/T 



Z = I DaDXexp { -y 



dr / d 2 r 



(32) 



[(Va) 2 + f(a 2 + .. + a%_,) + i\{o 2 - 1)} } , 

where a = {eri...<7/v} is the vector field, p s is the 
quantum-renormalized spin stiffness, / is the dimension- 
less anisotropy parameter (see main text, we consider 
here for simplicity the two-dimensional case a = 0). In 
the ordered phase (only this will be considered) the ex- 
citation spectrum in the zeroth order in 1/N, which is 
given by the poles of the unperturbed longitudinal and 
transverse Green's functions, contains a gap f 1 ' 2 for all 
the components a m except for m — N: 



G?(g) = (l 2 + f) , G°(q) = q- 



(33) 



The (dimensionless) self-energy of the transverse fluc- 
tuations is given by 



2T f d 2 q G?(k-q) 



N J (2n) 2 11(g) 



(34) 



where 



= I I 7§^°(p)G°(q + p) + Sg?( ? ), 



(27T) 
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and g = (N —1) / p s is the quantum-renormalized coupling 
constant. Note that Eq. (|3"4")l represents a contribution 



Hk/f 1/2 ) 



k 2 



N \\n(k/f 1 / 2 ) + 2ir<T 2 /(Tg) 



(37) 
)■ 



k 2 1 
— 0( - 

N V ln(fc// 1/2 



The first term in Eq. (|37j) can be transformed to the 
form, given by the renormalization-group approach (gen- 
eralization of the Eq. (IT71) to arbitrary JV); to this end we 
erform in Eq. ([37]) the standard replacement N — > JV— 2 
[l4^ . which is required because of considering first- 
order 1/JV terms only, and collect logarithm into the 
power. Collecting other terms to comply the lowest or- 
ders of perturbation result (fTTj) yields 



X + ~(q,0) 



JV- 1 
p s q 2 

X 1 + 

JV- I 
Psq 2 



^ln* 
2tt q 



l/(JV-2) 



Tg ln( 9 //V2) 



2tt 



a 2 + 



2tt 111 /T72 



^ln-2^ 



a 2 , Tg, _q_ 



-l/(JV-2) 



l/(N-2) 

(38) 



Replacement JV — > JV — 2 and collecting terms into power 
in the multiplier in the second line of Eq. Q38p is however 
unjustified near the magnetic transition temperature, 
since in that case this multiplier describes contribution 
of non-spin-wave degrees of freedom, which reduces the 
transverse susceptibility (N — 1)/JV times |14j. Instead, 
we find in that case, similarly to the two-dimensional 
case (3 



X + ~(q,0) 



N -I 



(39) 



JV- 1 



Y V2tt n / 1/2 , 



9 



Since the aboveconsidered calculation accounted for loga- 
rithmic terms only, the same results with the replacement 
/ — > r hold in the presence of the interlayer coupling. 
Comparing this to the spin-wave results © and ([TTT) one 



should have in mind that N = 3 in the physical case and 
the used zeroth-order 1 /N result for a 2 is analogous 
to the spin-wave result for S / Sq. 
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